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f^ ' The dynamic response of an interacting electron system is determined by an extension of the 

("^ ' relaxation-time approximation forced to obey local conservation laws for number, momentum and 

04 . energy. A consequence of these imposed constraints is that the local electron equilibrium distri- 

bution must have a space- and time-dependent chemical potential, drift velocity and temperature, 
flj ' Both quantum kinetic and semi-classical arguments are given, and we calculate and analyze the 

ry^ , corresponding analytical d-dimensional dielectric function. Dynamical correlation, arising from re- 

laxation effects, is shown to soften the plasmon dispersion of both two- and three-dimensional 
systems. Finally, we consider the consequences for a hydrodynamic theory of a d-dimensional in- 
teracting electron gas, and by incorporating the competition between relaxation and inertial effects 
we derive generalised hydrodynamic equations applicable to arbitrary frequencies. 
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The dynamic response of correlated and/or damped electron systems has both received and stimulated considerable theoretical 
work and continues to present challenges. The dynamic dielectric function (DF) of an undamped quantum electron gas was 
first calculated in the random phase approximation (RPA) by LindhardtJ. His approach assumed an infinite electron relaxation 
' . time, T -^ oo, and it was suggested that the effects of phenomenological damping could be incorporated by assuming a finite 
relaxation time. However it is well known that this single relaxation-time approximation fails to locally conserve the basic 
conservation laws for electron number, momentum and energy. This results in a number of incorrect experimental predictions 
such as alteration of the static DF even though damping is but a dynamic phenomenoih 

The first corrective step taken to rectify this situation was carried out by MerminB who was able to derive a DF which 
I conserved electron number during collisions. This was achieved by using a relaxation-time approximation in which the collisions 
^^ , relax the driven electron distribution not to its global equilibrium distribution, but to a local equilibrium distribution specified 
^ ■ by a local chemical potential, fi{r,t). This number-conserving approximatiisn has since been widely rused to study the effects 
00 ' of scatteringijn several different systems such as periodic crystal potentialsQ, small metallic particlesu and interacting storage 
\l , ring plasmasB. Nevertheless, the Mermin DF violates the two remaining conservation laws and, in a one-component system of 
C^ . electrons, any mechanism which induces relaxation of a non-equilibrium electron distribution may not violate any of the three 
C_^ ' laws. Even if, in the presence of external sources, momentum loss does occur then energy conservation may not necessarily be 
^— ' , affected as in the case of non-magnetic static impurities. In this paper we show how satisfaction of all three conservation laws, 
and various combinations thereof, may be achieved in determining the dynamic response of an electron gas. The proposed key 
to the solution is that the local equilibrium distribution must not only exhibit spatial and temporal variations in jj, but also in 
the drift velocity, v (to conserve momentum), and temperature, T (to conserve energy) ^-i This ideaO has, independently, also 
been recently introduced, and implemented in the context of generalised quantum liquidsLl. 

We begin by explicitly deriving the full conserving dielectric function (FCDF) as appropriate for a one-component, d- 

dimensional, quantum system of electrons. In this context the relaxation time approximation can be viewed as a way of 

^ , incorporating, at least approximately, dynamical correlation effects within the response of an jjiteracting electron gas. The 

i? ' corrections to the non-coUisional DF, arising from the conserving (relaxation-time) approximationu, thus serve the same purpose 

r^ ' as the much sought-after dynamical local field corrections of the DF of an interacting electron gas. Other straightforward 

applications of the conserving approximation, that will be discussed, include to the plasmon dispersion and conductivity of the 

electron system. 

1,^ We will also demonstrate that a correct accounting of relaxation effects is of considerable impei'tance in providing an accurate 

rS , hydrodynamic description of an electron gas. Though hydrodynamic models of an electron gasEI are popular by virtue of their 

S ■ simplicity, there remains a long-standing problem in that they may predict physically incorrect plasmon dispersion relationships 

in both two and three dimensions. The difficulties can be traced back to the assumption of adiabacity, valid only at low 

frequencies, where there is an overcompensation of relaxation effects. Inertial effects become important at higher frequencies 

and thus a stress tensor is required, in contrast to specification of the sole scalar pressure term used in traditional adiabatic 

hydrodynamics. This is particularly important in the case of three-dimensional plasmons where the strong Coulomb potential 

is responsible for a large frequency gap at zero wavevector. In contrast, relaxation effects from electron-electron scattering 

compete more strongly at lower frequencies to maintain local stress isotropy. It follows that a hydrodynamic model of an 

electron gasi generalised to arbitrary frequency, must successfully incorporate competition between both effects. Tokatly and 

Fankratovtj have also recently succeeded in deriving such a model for a three-dimensional electron gas based on the Boltzmann 

transport equation, and have gone on to show how such generalised hydrodynamics can be applied to the Landau theory of 

Fermi liquids. Here we present a derivation of the model for a d-dimensional electron gas starting from the microscopic quantum 

dynamics but with the relaxation-time approximation forced to conserve all three laws, and by way of example, demonstrate 

the manner in which plasmon dispersions are then correctly deduced. 
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The plan of this paper is as follows : in Sec. nj we introduce a relaxation-time approximation within a formalism of quanturfl 
transport theory which allows for local chemical potential, drift velocity and temperature variations. In Sec. |l| we develop 
applications of the conserving relaxation-time approximation by, firstly, deriving the new DF and comparing its properties, 
including plasmon dispersion, with previous models, and secondly, considering the conductivity. In Sec. [[V| we show how an 
analogous semi-classical derivation can be made, and in Sec. M we derive a hydrodynamic model of an electron gas generalised 
to arbitrary frequencies. 

II. QUANTUM KINETICS 
A. Model 

The model we consider, primarily, is a system of charged electrons embedded in a neutralizing background (jellium model), 
whereby electrons are only scattered by other electrons, and the dynamics of such scattering events are constrained by all 
the conservation laws. In any realistic setting there will also be phonons and impurities giving rise to additional scattering 
mechanisms with fewer conservation constraints, and indeed the loss of momentum and energy is important in a discussion of 
the conductivity and dissipation of the system under the influence of an external field. If all possible scattering processes are 
treated at the same phenomenological level we need only to discuss in detail the more constrained problem of a one-component 
system (i.e. electron-electron scattering only) to understand how conservation constraints may be applied in general to the 
other scattering mechanisms. The one-component model has the additional virtue of allowing us to calculate dynamical local 
field corrections of the DF arising entirely from electron-electron correlation effects. 

The treatment of a d-dimensional many-electron system is simplified by reduction to an underlying single particle description. 
In practice this means tracing out all the degrees of freedom associated with the other remaining electrons and embedding the 
many-body interaction term in the dynamics of a single electron. In this way, the dynamics of a single electron is governed by 
two terms : i) an effective single particle Hamiltonian, and ii) a term describing interactions with the other particles. 

B. Relaxation-Time Approximation 

Hence, we begin by considering the dynamics of the quantum one-electron statistical operator operator p. The fermionic 
single electron equilibrium statistical operator is given by {h = 1) 

po = ^ ^ ^ , (1) 

g{Ho-k'V-M)/fcBT _|_ I 

where Ho and k are the Hamiltonian and the momentum operator for an independent and non-interacting electron. Henceforth 
we shall refer to po as the global equilibrium statistical operator. The distribution of a moving gas differs from that at rest by 
a Galilean transformation of the drift velocity v = (k)/rra, which is assumed to be much smaller than the Fermi velocity. With 
a view to future simplification we set the global equilibrium drift velocity, v^,, equal to zero (i.e. the change from equilibrium 
and the magnitude of the drift velocity is due to an applied field). Note also that in a crystalline lattice the mass of an electron, 
m, which enters the calculations is taken to be the effective mass thereby incorporating, at least approximately, band structure. 
For any eigenstate |i) of Ho, 

where fi is the Fermi-Dirac occupation function for that state. Upon application of an external perturbation, C/oxt, and ensuing 
internal response, Uint, the dynamics of the total single particle statistical operator, p, is governed by the quantum Liouville 
equation 



a = ft) . (3) 



f + -i«— u,; 



' coll 



where H is now the effective Hamiltonian, i.e. 



H = Ho + U, (4) 

U = f/cxt + Ulnt, 

and the term on the right of Eq. (bl) is the collision term. Note that H is not the full Hamiltonian since it represents mean-field 
(Hartree-like) dynamics containing only one-body operators, whereas the two-body terms are subsumed in the collision term. 
The standard relaxation-time approximation is now invoked, namely 

'dp\ _ p - pio, 



fdp\ _ p - pioc . , 



The usual interpretation of the approximation is as follows : the electron ensemble evolves in time as determined by the 
left-hand side of the Liouville equation, Eq. (bl), but in a time interval 5t a fraction, St/r, of them collide and are redistributed 
according to a local equilibrium statistical operator, pioc- The energy dependence of r can be ignored if we assume that the 
electron levels involved in the collision processes lie close to the Fermi energy. 

However, the familiar approximation of setting pioc equal to po leads to a violation of conservation laws. To determine the 
correct (dynamical) local equilibrium state we observe from classical thermodynamics that any state of an ensemble of particles 
can be specified by five parameters (in three dimensions), two coming from the thermodynamic variables /x and T, and three 
from V. These parameters are well defined quantities (constants) in a global equilibrium state, as specified by po- If we now 
allow these thermodynamic parameters to become functions of space and time, such that there is no entropy production, then 
the entire electron gas is said to be near-equilibrium and locally the electron distribution is given by the form of the global 
equilibrium distribution, po, but now with local parameters. These parameters must be constrained, and hence determined by, 
conservation laws for average particle number, momentum and energy. In other words, by allowing local fluctuations of p, v 
and T we gain the necessary degrees of freedom to satisfy the conservation laws. The thermodynamics of the local equilibrium 
state and the physical meaning of thermodynamic quantities in a near-equilibrium gas is further elaborated in Appendix N. 

As appropriate for linear response theory we treat the varaiations of the local equilibrium state as small, and thus the local 
Hamiltonian speciflying the local distribution, can be expanded as 



pioc(Hioc) = ?o{Hq) + Api + O(A^), 



^loc ~ Ho + A 






+ 0(A^ 



(6) 
(7) 



where A is linked to the order of the expansion. By definition, the local equilibrium density operator must commute with the 
local Hamiltonian and thus we obtain the following matrix element to first order (after setting A = 1), 



(k|pioclk') = /k(5kk' 






5p(k - k') + ^ . 5v(k - k') + fil±|^^5T(k - k') 



(8) 



Note that we have implicitly assumed the electron state distribution to be isotropic and homogeneous in k— space. To determine 
the response we assume that the perturbing potential, and hence other dynamic quantities, varies as e*'''''^"*\ We restrict 
our analysis to a longitudinal response such that Vf/cxt x q = 0. The dynamics of the first order density matrix, 5p, is then 
determined by Eq. (fBT) and Eq. (H) giving 



a;(k|<rplk') - (ek - ek')(kl<rplk') + (f^, - /k)f/(k - k') 



(9) 



T 1 Ek - Ek' 



5p(k - k') + ^ . 5v(k - k') + £il±^^5T(k - k') 



It is convenient to change variables, k -^ p -|- iq and k' ^ p — iq, thereby making apparent the disturbance wavevector, q. 



(a;-A6) (p-fiq|5>|p-iq)+A/C/(q). 
-- (p+-qlMp--q) + ^ 



where Ae = 
respectively. 



^Pfiq - Viq ^nd A/ = /, 



p^fq 



-/. 



(10) 



Mq) + P-Mq)+(^ + Viq-^^'^^''^ 



2 P-2I '^ J T 

are the single-particle changes in the energy and distribution function 



Conservation Laws 



The expectation values for the conserved quantities (number, momentum and energy), in a particular distribution, are 
obtained by taking the appropriate moments of the distribution in momentum-space. In general, momentum and energy will 
not be fully conserved in the presence of other degrees of freedom, e.g. phonons and magn etic im purities. In fact, loss of 
momentum is necessary to render a finite conductivity and this will be discussed later in Sec. [II C Here we will concentrate 



on a full-conserving model and to ensure that all these quantities are collision invariants we must obtain the same expectation 
values before and after each collision. The formal representation of this statement is 



(27r)« 



c( P (P + 2'^\P^ ?1°<;|P - 2^) 




(11) 



Note that the quantities being conserved by Eq. (Illl) are statistically-averaged, i.e. it is the expectation values of the density, 
momentum-density and energy-density which are being conserved. Of course, angular momentum is also conserved, but in 



a non-rotating system of point particles we will not obtain an additional constraint from this conservation law. This can be 
seen from the fact that the angular momentum operator is merely derived from a bilinear combination of the position and 
momentum operators, and thus conservation of angular momentum is automatically enforced. 

It follows that the taking of moments of either side of the linearised equation of motion, Eq. (|lC| ) , must necessarily give zero. 
We now proceed to describe the dynamical response of an interacting electron gas where the free part has an energy dispersion 
given by i.e. -ffolp) ~ ^pIp) = P^/27n|p). Enforcement of the conservation constraints on the right-hand side of Eq. nUM gives 



where :r is a column vector of the small changes in number density, momentum density and energy density, i.e. 



Further, 



and 



M = 



gj" (2^/'^'P<P+^'ll^'^^IP-^^>(p j 



-Bo {Bon-B2/2m~Boq^/8m)/T 

no 

-B2/2m (B2M/2m-B4/4m^-B2gVl6m^)/r 



/ 




(12) 



(13) 



(14) 



(15) 



For brevity we have also introduced a function, B„, the n^^ momentum moment of the integrand of the static Lindhard 
(free-electron) polarizability function, namely 



Sn(q) 



(2TTy 



dp\P\ -^, 



It will also be useful, in the ensuing analysis, to introduce the following related dynamic functions, namely 



C„{ci,uj) = 



dp\p\' 



A/ 



and 



D„{ci,Lu) 



{2tiY J ''"'' Cj-Ae 

iCn — UJrBn 



(16) 



(17) 



(18) 

uj I 

where Cj — uj + i/r. Analytic expressions for moments of the dynamic d-dimensional Lindhard polarizability function in the 
limit |q| ^ A:f are evaluated in Appendix Fs. 

To solve for / we must obtain a second set of equations, obeyed by x, Eq. (|l3), by taking the relevant moments of Eq. (Mm. 
We can proceed in either of two related ways. The first, is to solve for {p+^ci\5p\p~^q) in Eq. (lid) and then take moments, but 
this results in rather long and cumbersome intermediate expressions. We outline such a generalized formalism in the Appendix 
O Second, as will be seen below, we can partially ease the task by noting that the eventual outcome of all these calculations 
is to obtain a self-consistent expression for 5n. In view of this we need only be concerned with the task of finding alternative 
expressions for the last two components of x. By imposing conservation of number on the left-hand side of Eq. mjf) we obtain 
an equation of continuity, from which we have 



5i = Sn^-q^. 



(19) 



However there exists no similar strategy for obtaining a simple expression for S£ as a function of 5n and we must resort to the 
more formal first method, as outlined in Appendix O This leads to 



S£ 



-Uici, 



Ci , B2 + C2 



2m 



2miijjT 



5 II 



C2 



2mq'^iT 



— 5kq + 



B4 + C2 



4m2T 



-f 



g 

8m 



M 



B2 + C2 



2mT 



ST. 



(20) 



After eliminating the elements of x the local thermodynamic variations are then found to be 

Boq^\ ( icjC2 



5/i(q,w) = 



Bo2mfi — B2 



rq^no 



5n + UC2\ +Sn( 022771^-^^-04] 



-D4-B0 — O2B2 

(5v(q,a;) = <57i— ^q. 



and 



(21) 
(22) 

O4B0 — O2B2 

These are among the primary results of this paper; they have immediate application in the determination of the dynamic DF 
for a relaxing (correlated) electron system. 



ST{ci,uj) = 2mT- 



5n[02 + 



7ioq- 
iujBoC2 



UC2B0 



III. APPLICATIONS 



Dielectric Function 



By substituting the ex pre ssions for the thermodynamic fluctuations, Eq. (Ell), Eq. ( |22[ ) and Eq. (|23[), in the linearized 
transport equation, Eq. ([lO|), we can solve for the density perturbation, 5n. To obtain a self-consistent expression for the 
response we specify a Hartree potential for the internal repsonse, i.e. f/int(q) ~ —5n{q)V{q) where ^(q) is the Fourier- 
transformed Coulomb potential, 



v{ci) = 2'*-V('*-^>/^r[i(d-i)]eV 



(24) 



^(q) 



47reVg^ 
27ve^/q, 



(d=3), 
(d=2). 



The longitudinal dielectric function follows immediately from linear response theory, namely 

<57i(q, Lu) 



e{ci,cu) = l + V{ci) 



c/(q) 



The expression for the fully conserving dielectric function (FCDF) is then found to be 



6=(q,c^) = l + y(q)- 



where 



and 



1 + F ' 
C2 \ C2B0 — C0B2 



F = 



Coir J D4B0 — D2B2 

D2C2 — -D4C0 H 5 — (C2-B0 — Co-B: 

rq^'no 

-04-80 - D2B2 



- 1 



+ 



iujCo 



(25) 

(26) 
(27) 

(28) 



are the conserving damping corrections. Though cumbersome, some important properties of the degenerate electron gas are 
fairly easy to deduce from Eq. (|26[), Eq. (ETI) and Eq. (GSl) with the aid of analytical expressions for certain integrals (Appendix 
g). First, by taking the limit r ^ 00 in the FCDF we reassuringly recover the celebrated Lindhard dielectric function in the 
absence of damping, i.e. 

2 /■ .. fu 



6"(q,a;) = l + V(q)- 



d>- 



'p+-q 



-/p-iq 



(29) 



' {2nY J "'a;- p- q/m + iO+' 

Second, by taking the static limit a; ^ we find 

6=(q,0) = e°(q,0), (30) 

as expected since relaxation processes ought to have no bearing on the static properties of the electron gas. Hence we also 
obtain the correct Thomas-Fermi screening result at low q, i.e. 



gTF 



^2d-2^(d-i)/2p((^_^)/2)noe 



i/{d-i) 



(31) 



where Ep ~ m ^2^ ^'''7r[r(l + rf/2)]^'''?i^''' is the d-dimensional Fermi energy. Third, in the long wavelength limit, q ^ 0, we 
retrieve the same form for e'^ as that obtained in the absence of damping, namely 



lim e''(q,a;) = 1 1, 



q^O 



where 



2 noq 



^'r. = —V{<\), 



(32) 



(33) 



as required by the sum rules for an electron liquidllil. 

Using the general form of the dielectric function as given by Eq. (|2^) we may immediately make comparisons with other 
prominent models proposed in the literature : (a) the undamped Lindhard DF (RPA) corresponding to the choices E = Q, F = Q 
and r -> 00, (b) the damped Lindhard DF: -E = 0, F = and finite r, (c) the Mermin DF : _B = 0, F = [i/(w-r)][l -|- Co/Bo] 
and finite r, and finally (d) the FCDF given by Eq. (ETI), Eq. ( |2q ) and finite r. The real and imaginary parts of the various 
dielectric functions are shown for comparison in three dimensions (Fig. hi and Fig. bl) and two dimensions (Fig. H and Fig. M). 
It is not surprising that in all cases the behavior of the FCDF most closely resembles the RPA as it is only in these models 
that all the conservation laws are enforced. 



B. Plasmon Dispersion 

A simple application of any proposed DF is to determine the dispersion relationship, ct;(g), of bulk longitudinal plasmons 
by imposing the condition : e(g, uj) = 0. Plasmon dispersion data can be obtained by scattering experiments involving X-rays 
or electrons. The measured quantity in these experiments is the energy loss of the scatterer which is directly related to the 
imaginary part of the DF. For the FCDF in the long wavelength limit it can be shown that 

lim Im(e^"'") — —ujpS{uj — Wp). (34) 

q—*0 2 

Herfe« plasmons are the only energy loss at q = and we have exact fulfillment of the long wavelength "perfect screening" sum 
ruleO: 

f°° 1 _i n 

lim / — Ime; {q,uj)dLO=—. (35) 

'^"Jg CO 2 

We analyze the solutions for plasmon dispersions arising from the various dielectric models as listed in the previous section. 
For three dimensions the differing models are compared in Fig. H. As expected, the softening of the plasmon frequency from 
both the FCDF (at finite q) and the Mermin DF increases with r^ and 1/r. The salient differences between the models are 
exhibited in the limit q —> 0. The FCDF gives rise-ito the same plasmon energy gap, ujp, as the undamped Lindhard DF, as 
required from the sum rules for an electron liquidllj. In Ref. H the Mermin DF was apparently deduced to obey the perfect 
screening sum rule. However this claim can only be true in the limit w — » cxa (or r — > oo) and it must fail at finite lo as the 3D 
plasmon dispersion indicates. In fact, after a little algebra, it is possible to show that the Mermin DF gives rise to a Lorentzian 
spread in the energy loss, 

2 ; 

limlm(eMormin) = -7-^ 2^1/ I ^2 ■ (36) 

The plasmon dispersion for two-dimensional systems is plotted in Fig. g where it should be noted that t is assumed to be 
large enough such that localisation effects are negligible. As in the three-dimensional case, the plasma excitations from the 
Mermin DF and the FCDF are slightly softened with respect to the plasmon from RPA. Interestingly, the Mermin DF predicts 
a wavevector culrfctpff below which a plasma mode cannot exist. This cutoff, which increases with r^ and 1/r, has also been 
previously notedta via the summation of ladder diagrams in the disorderd polarazibility function. The q — > limit of the 
FCDF matches that of the RPA DF in agreement with the aforementioned sum rule. For comparable values of rs and r the 
FCDF shows a greater deviation away from the RPA DF in two dimensions than in three. Clearly, scattering processes are less 
important for the three-dimensional plasmon because of the large frequency gap. 

The damping of the two-dimensional plasmon, either from the presence of external impurities (Mermin DF) or electron- 
electron scattering (FCDF), occurs before the onset of Landau damping (not shown in figure) and thus may be amenable 
to experimental verification. In fact it has already been appreciated that experimental studies of two-dimensional electron 
systems are an important test of theoreticait4pflatments of many-body effects (beyond RPA) which have been shown to be more 
important in two dimensions than in threecStJ. 

C. Conductivity 

In a one-component system, and in an approximation where momentum is rigorously conserved, we expect a steady state 
D.C. current to increase indefinitely upon application of a constant external field. It can easily be shown that the conductivity, 
o, as determined by the FCDF, is indeed infinite for the full-conserving model. To accommodate loss of momentum, thus 
rendering a finite conductivity, we return to the equation of motion 

(p+oql'5p|p-oq)('^-p-q/"T-) + A/c/(q) = — (p+T^qlp-piodp-Trq). (37) 

II T I I 

and impose static (a; — > 0) and homogeneous (q -^ 0) conditions, giving (in real space), 

5i = <5jioc - norVC/, (38) 

where 

<5Jioc = j—y j dpe*<''"'""'^'(p+-qlpioc - P()lp--q)p. (39) 

If a fraction. A, of the momentum density is conserved during collisions, i.e. 

<5jioc = A5j, (40) 

then, using e'^Jj/m = — ctVC/, we obtain the familiar Drude conductivity 

a = noe Trn/m, (41) 



where we have defined an effective momentum-relaxation time, 

rm = r/(l-A). (42) 

As expected, the conductivity diverges in the conserving hmit, A -^ 1. In the presence of external scatterers (e.g. fixed 
impurities or Umklapp scattering), the other limit, A — > 0, corresponds to isotropic scattering whereby the net momentum 
density is totally destroyed during collisions, and intermediate values of A refiect the degree of anisotropic scattering. 

We may also include a second contribution to the collision integral stemming from additional scattering processes. These 
may include inelastic processes which we have up until now ignored. In other words, ther^levant conductivity relaxation time, 
Tcond, niay contain a contribution from both r and an often larger energy relaxation timetj, r'. To demonstrate this we assume, 
but without any strict justification, a (decoupled) relaxation time approximation for the inelastic contribution to the collision 
integral, 

\dt ) coll r t' 

where pj^c describes a state of local equilibrium where all energy and momentum from the motion of particles has been lost. 
The expression for the conductivity then becomes. 



■ noe'Tcond/m, (44) 



where t^J^^ = ''"m ""^ + ''"'/ ^ (i-e. Matthiessen's rule). 



IV. SEMI-CLASSICAL LIMIT 

We will now derive the dielectric function in a fully classical formalism except that at the end of the calculations we will still 
be free to employ a distribution appropriate for either a classical or a fermionic gas. The emerging result for the degenerate 
fermionic electron gas, whose equilibrium distribution is governed by quantum statistics, will be what we call the semi-classical 
limit. 

The central dynamic quantity in classical transport theory is the phase space electron distribution /(p, r, t) governed by the 
Boltzmann equation. In the presence of an electric field, E, the equation of motion in the relaxation time approximation is 
given by 

§f+P.§/+eE.^ = -i^^. (45) 

ot m or op T 

Magnetic effects are of second order and so disappear when the equation of motion is linearized. The total electric force, eE = 
—VU, has a contribution from an external source and an internal contribution from the induced charge density perturbation, 
as discussed in the quantum case. To deternune the time development we now assume that the distribution function is slightly 
perturbed from equilibrium, and write 

/(p,r,i) = /o(p) + 5/(p,r,i). (46) 

As in the quantum case we work within the grand canonical ensemble which serves to define the global macroscopic temperature 
and chemical potential. To satisfy the local conservation laws, the local equilibrium distribution must be different from the 
global equilibrium distribution, /o, and to first order is given by 

where 

dpP/. (48) 

m 

Again, in our model we assume that the electrons have no initial drift velocity and only experience a net drift upon application 
of an electric field. Rewriting the local distribution as 

fioc = fo + -Q^5fi+—p-5v+— ^ ST, (49) 

makes it transparent that it is indeed just the diagonal element of the corresponding quantum statistical operator in the limit 
q — > 0. The Boltzmann transport equation, Eq. (f45|), now assumes the following form when linearized. 



dt m dr dp t t de 



Sf, + p-Sv+^-^j;^ST 



(50) 



As in Sec. II C we now examine the consequences of applying the conservation laws. The terms 5fi, Sv, and ST are independent 
of the electron momentum and must be determined by the local conservation laws for the collision invariants : average number, 
momentum and energy, i.e. 



d%{f - /loc) \ P i = < I . (51) 

[ p^/2m J i J 

Proceeding in a manner similar to that described in section O, the solutions are, 

(B'o2mfi - B2) I ^^^(Sn + UC2 ] + Sn {D^^'imfi - D'^) 



Mq,^) = ^ n, n, n, n, ' (52) 



(5v(q,tj) = Sn -q, (53) 

and 

Sn (d'^ + 'J^^\ + UC^B'o 

We have introduced the classical versions of the integral functions defined in section [11, i.e. 

^"-(2^/^^^iPl-f' (^^) 

(ztt)" J dp Lj — p ■ q/m 



and 



, _ C'„ + iuJTB'n 



ILOT 



(57) 



Note that the expressions for the thermodynamic variations also differ from their quantum counterparts in that terms of order 
(f are absent. 

We can now obtain the dielectric response of the system with a treatment similar to that given above. The DF that follows 
from solving for Sn in the transport equation is identical to Eq. (Gq) except that the integral functions, i3„, Cn and D„ 
are replaced by their classical counterparts, B'„, C'„ and D'„. The interesting point here is that functional form is therefore 
preserved in both the quantum result and in the classical approximation. Interestingly, an early attempt to "guess" a quantum 
number-conserving response from the classical response led to erroneous resultsEJ. However the work presented in this paper 
suggests that the the correct response may be correctly guessed from the classical response but only if the classical response 
can be written in terms of macroscopic moments of the Lindhard polarizability function. 

The semi-classical result follows by taking /o to be the equilibrium fermionic distribution. In this case the integral functions, 
B'n, CL a nd D'„, a nd he nce the semi-classical DF, represent the low q expansion of their quantum (unprimed) counterparts in 



Sec. [IC and Sec. Ill A, respectively. 



V. HYDRODYNAMICS OF A DEGENERATE ELECTRON GAS 

A hydrodynamic description is obtained by coarse-graining over the microscopic degrees of freedom in such a way that the 
dynamical equations are expressed in terms of macroscopic quantities only. The hydrodynamic equations are traditionaHji 
derived from the classical Boltzmann equation Eq. ( |45| ) and this indeed is the approach employed by Tokatly and Pankratovlij 
who also recently formulated the hydrodynamic theory of an electron gas to arbitrary frequencies. Here we provide a derivation 
from the quantum description of the microscopic dynamics and begin by writing the macroscopic quantities in terms of the 
one-particle statistical operator, 

n(r,i) = (rl^r), (58) 

j(r,t) =mn(r,t)v(r,t) = i{rj{p, p}lr), (59) 

U{r,t) = {r\U\r), (60) 

and 

-P»j(r,i) = 7^— (r!{p, {Pi - mvi){pi - mvi)}\r), (61) 

Am 

where {, } denotes symmetrization (i.e. anticommutation) , n and j can be shown to be equivalent to the definitions in Eq. ([13|), 

and Pij is the local stress tensor in the comoving frame. Higher momentum moments can be neglected in a second order 

theory of a degenerate electron gas, i.e. thermal effects are negligible since ksT <^ Ef- The conserving constraints can also be 

expressed in a similar manner, namely 
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{r|p-pioe!r)=0, (62) 

(rl{(p-pio,),p}!r) = 0, (63) 
and 

(r|{(p-p,o,),i?o}|r)=0. (64) 

The hydrodynamic equations follow by taking diagonal elements (in position representation) of the quantum dynamical 
equation, which we repeat here, 



dp . 



2m 



^i^^^, (65) 



but symmetrized with the relevant operators from each conservation law. Before proceeding we note than an alternative 
approach to deriving the equations of motion of thermodynamic coordinates from the quantum LiouviWq equation, but not 
necessarily close to equilibrium, has been discussed via the time-dependent projection operator techniquell3. 

The first hydrodynamic equation, the equation of continuity, is simply obtained by taking the diagonal element of Eq. (p5[), 
eventually giving 

Dn 

— + nV^v, = 0, (66) 

where Dt = dt + I'iVi is the covariant derivative. The derivation requires use of the fact that for a one-body potential, 
(r| [17. ci\\r) — and that {r| [p, A]\r) — i\7{r\A\r) for any operator A. The next hydrodynamic equation follows by anticommuting 
Eq. (p5|) with p and then taking the diagonal element, culminating in 

mn^ + ViP + Vjni-j-nV.U = 0, (67) 

where the symmetric stress tensor, Pij , has been decomposed into scalar and traceless parts, i.e. 

Pij = pSij + TTij (68) 



and 



The scalar d-dimensional pressure of the gas in local equilibrium is given by 



(69) 



P = ^TrP.^ (70) 



which for an ideal and degenerate Fermi gas, of spin degeneracy 2, can be calculated to be 



fe.T^^^^2--/Mr(l + W^-nr^^^//A^X-. (,,^ 



V d+2 m \\ Ef 

where Z is the grand canonical partition function. The third hydrodynamic equation, obtained by anticommuting Eq. ( pq ) 
with piPj , is given to second order by 

^rr + P^l^kVk + Prk^kV, + Pk,\7klH = -^. (72) 

Ut T 

By decomposing Pij into p, Eq. ([7l|), and -Kij we actually get two equations from Eq. ([7^), namely 

TTik^kVi = (73) 

and 

— — ^ +nij\/kVk +'!^ik^kVj +7^kj^kVi — -i'5ij7r,ifcVfcU^ +p IViVj + VjHi — -SijVkVkj = —■ (74) 

The complete set of hydrodynamical equations is then specified by Eq. (pw, Eq. (pTI), Eq. (0), Eq. ( [73| ) and Eq. ([74|). 
Traditional Bloch (adiabatic) hydrodynamics is recovered by setting -Kij = 0. The presence of nij stems from the anisotropy 
of the stress tensor, i.e. the inability of the dynamic electron gas to achieve local equilibrium. The content of Eq. (p4J) shows 
that the extent of anisotropic dynamics are tempered by relaxation effects which drive the system back to an isotropic state of 
local equilibrium. 

We now demonstrate the efficacy of these equations in correcting the shortcomings of the traditional adiabatic hydrodynamics 
in predicting the longitudinal plasmon dispersion. The simplest route to calculating the dispersion of the collective mode, in 
say the x direction, is to first linearize the hydrodynamic equations, giving 

dSn „ „ ,„^, 

-^ + noVxVx = 0, (75) 



dvx 
mno^— +VxSp + Vx5Trxx — noVxU = 0, 
ot 

d no 



and 
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(76) 

(77) 

and then to consider a plane- wave solution e''-'^"'^*' in the absence of an external field, i.e. U = Uint- Ignoring retardation 
effects, the potential energy U{r, t) is simply the convolution of the density perturbation and the Coulomb potential energy, 
which can be expressed in Fourier space as U{q) = —5n{q)V{q), i.e. mean-field Hartree potential. Solving the linearized 
equations then gives the following solution, 



dS-K, 



OCX ^(-L ^ 

H ; Po^/xVx = 



S-Kx 



UJ^ = uliq) + I3\u)q^ 



where ijJp{q) is given by Eq. (BSf) and 



r(^) 



I3{uj) 



rd + 2 2d-2 /cj 
d d \uj 


)] 


Po 
mno 


'l 2d-2 (u\ 
d^ d{d + 2) V(D/ 


1 


/2 

VF- 



(79) 
(80) 

(81) 



The coefficient /3 interpolates between the correct limiting forms of both the hydrodynamic regime {lot ^ 1), ini-agreement 
with conventional adiabatic hydrodynamic theory, and the coUisionless regime (tor ^ 1), in agreement with RPAuZI. In a d- 
dimensional electron gas, the low frequency limit entails the fact that collisions are effective in maintaining a local equilibrium 
and thus a conduction electron possesses d translational degrees of freedom. In the opposite high-frequency limit the influence 
of collisions is negligible and the electron has effectively one translational degree of freedom, corresponding to motion in the 
direction of the applied field. The notion of effective number of degrees of freedom, /, can be made quantitative by comparing 
Eq. (teol) with the same result in the adiabatic limit, i.e. ignoring the second term in brackets, giving 



The limiting forms of /(a;) and I3{uj) may be summarized as 



du} 



LU + {d — l)tJ 






{UJT < 1) 



(82) 



(83) 



and 



(lot > 1) 



(84) 



fH = 1, 

l3{Lu) = [3/{d + 2)]'^\p. 

In three dimensions, Eq. (pll) and Eq. ( p2[ ) agree with those derived by HalevitJ who determined the u-dependence of f3 by 
comparing a (damped) hydrodynamic DF with the Mermin DF. The need for ar^equency-dependent j3 has also been noted in 
the context of a hydrodynamic description of inhomogeneous electron sy*fe«mai3. The cross-over between the coUisional and 
coUisionless regime has also recently been discussed in an interesting papeiE3 discussing the viscoelastic behavior of an electron 
liquid. 

Finally, we noter-fehat a popular way of implementing phenomenological damping in hydrodynamic models of an electron gas 
through the use ofEj 



Dn 



+ nViVi = 0, 



and 



Dvi ViP 



V,U- 



(85) 



(86) 



Dt n T 

is simply the hydrodynamic approximation of the Mermin model in the adiabatic regime and where all momentum is destroyed 
during collisions, i.e. Eq. (p5|) and Eq. (feq) are obtained from Eq. (p5|) with the following constraints, namely 



and 



{r|p-pioclr) = 0, 
TTij = 0, 

{rl{pioc,p}lr)=0. 



(87) 
(88) 

(89) 



VI. CONCLUSION 11 

The relaxation-time approximation has been applied to the quantum dynamics of a rf-dimensional electron gas where (average) 
number, momentum and energy are conserved during collisions. It has been shown that the requirement of conservation of 
all collision invariants necessitates a state of local equilibrium, to which the electrons relax, with a space- and time-dependent 
chemical potential, drift velocity and temperature. The ensuing dielectric response has been determined and compared with 
others in the literature, revealing that, in general, that imposition of the conservation laws tends to make the dielectric response 
in the relaxation-time approximation more akin to the RPA dielectric response than without them. The FCDF attains the 
correct static limit and obeys the perfect screening sum rule. 

The resultant plasmon dispersion has been compared with others and relaxation has been shown to induce softening of the 
plasmon with increasing rs and 1/r. In both three and two dimensions and in the long wavelength limit, the plasmon dispersion 
from FCDF is in good agreement with that derived from the RPA. In general, correlation effects are more pronounced in two 
dimensions than in three and thus further experimental study of two-dimensional plasmons ought to provide a significant test 
of models of dynamical correlation. 

The derivation of the conserving reponse has been repeated but in a semi-classical formalism using the Boltzmann transport 
equation. The functional form of the ensuing response is then identical to that obtained by quantum kinetic arguments. 

The hydrodynamic equations of a degenerate electron gas are derived from the quantum microscopic equation of motion 
using the relaxation time approximation and fulfillment of all the conservation laws. The departure of the stress tensor from 
a scalar is due to inertial effects but counteracted by relaxation effects. The competition between these two effects gives rise 
to a generalization of the hydrodynamics applicable to all frequencies and thus correctly reproduces the collision-dominated 
adiabatic limit [ur <^ 1) and the coUisionless limit {ur ^ 1). 

Given that the FCDF aims to incorporate dynamic electron correlation effects via a phenomenological r it is interesting to 
compare this with dielectr jei funtions which include correlation effects via a many-body field-theoretic approach giving rise to 
dynamic local field factor^. In principle the two different approaches ought to exhibit the same qualitative behavior (with 
respect to differences from the RPA) and indeed this is what we observe. In this context we can view the relaxation-time 
approximation as being a crude way of including Feynman diagrams beyond the RPA in the proper polarizability of an electron 
liquid. To extend the comparisons in a more quantitative and direct manner further work is required, i.e. calculation of an 
expression for r and inclusion of exchange in the FCDF. 

Finally, we repeat here that the analysis in this paper explicitly calculates the response function only when all the conservation 
laws are obeyed, as required for a one-component system, but, using the results of Appendix O, the response functions can be 
calculated which obey fewer and in different combinations. This is the case required in, for example, two-component systems 
or when crystalline effects are incorporated by way of Umklapp scattering. 
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APPENDIX A: THERMODYNAMICS OF A LOCAL EQUILIBRIUM STATE 

In this appendix we elaborate on what is meant by local values of ^, v and T, in the thermodynamic sense. 

We first recall how temperature and chemical potential are defined for a large (macroscopic) system. A, lying in contact with 
a heat reservoir, B, with which it can exchange particles and energy. Thermodynamic equilibrium of the system-|-reservior 
setup is achieved when the entire entropy of the system-|-reservoir is maximized. 

From the conservation of total particle number, energy and volume we arrive at the equality of three thermodynamic variables 
which we identify as the temperature {Ta = Tb) and the chemical potential (^a = A*s)- Throughout the system ^ and T 
are constants independent of space and time. By appealing to the strong Coulomb attraction between the electrons and the 
positive background we can infer that the velocity field, va, and hence the pressure, is also constant throughout. The constant 
nature of the thermodynamic variables is due to the effect of relaxation processes which drive the system towards equilibrium. 
Now, if the system is perturbed slightly, these variables may attain a spatial variation and will still remain thermodynamically 
well-defined if the subsection of the system over which they vary slightly is much greater than the average distance between 
collisions. Thus we require that 

Iq^mt < 1, (A2) 

where q is the wavevector of the spatial fluctuations of ^, v and T, and /mf is the mean free path which can be approximated 
as vpT. In addition the time variation of the thermodynamic variables must be slower than the relaxation time scale in order 
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for the relaxation processes to maintain the state of local equilibrium, i.e. 

ujT < 1. (A3) 

This constraint is consistent with the fact that the scattering term, and hence the presence of /loc only becomes important in 
the limit that r is small. 

These considerations suggest that the local equilibrium state may be an adiabatic one and indeed we can show that /loc is a 
small perturbation away from /o such that there is no entropy production. We begin by recalling the Von Neumann definition 
of entropy, 

S = -fcBTr(pgclnpgc), (A4) 

where the trace is carried out over all possible microstates in the grand canonical ensemble and pgc is the statistical operatpj 
of the grand canonical ensemble. This expression can be written in terms of the one-particle fermionic statistical operator, pp3 

S = -fcBtr[plnp+(l-p)ln(l-p)], (A5) 

where tr denotes the trace over one-particle states. In a non-equilibrium distribution the occupation and distribution of states 
acquire both a spatial and time dependence and thus, in general, the entropy will be a function of time. The time variation of 
pioc is determined by substitution into Eq. ( pq ) giving 

^=-i[i?,pioe], (A6) 

and it follows that when p = pioc 

We note that this derivation is essentially the quantum analogue of Boltzmann's Q (or 7i)-theorem. Thus the variations of 
/i, V and T, that characterize the local equilibrium state, are constrained to prevent the production of entropy in the system. 
The situation is identical to the propagation of sound in an Euler fluid or flrst sound in superfluid Helium II (to lowest order in 
{Sn/ST)p): a linearized disturbance away from equilibrium gives rise to an isentropic (constant entropy density) longitudinal 
traveling wave with local variations in density, pressure and temperature. 

APPENDIX B: USEFUL INTEGRALS 

The integrals defined in this paper can be expressed as moments of the polarizability function, x, in the Lindhard dielectric 
constant. The Lindhard polarizability, in d dimensions, is given by 

(271)" J Lo — p ■ q/m + «0+ 

We now define the n— th momentum moment of the polarizability function as 



2 / jd I i?i •'p+iq ■'p-fq 



(2-k)'^ J cj — p • q/m -I- iO+ 

Here we demonstrate how analytical expressions may be obtained in the low q limit for a degenerate electron gas. In this 
limit, 

Rex4q,^) = jAl^ [d'p^ '"'"^'"^ , (B3) 

(zttY^ J ae ujm — p • q 






Imx(q,tu) = ■;^— ^j^ / dp^lpl"p-q5(a;m-p- q). (B4) 



In spherical coordinates the real part becomes. 



Rexn(q,t^) = ~~T2§)r'^ j '^P'^P'^ ^^" / d(t> \u)I[(t>,p] ~ —sm^ ^0j , (B5) 



m,p\ = 



cos4> 



where Q,d = 27r'^"/r(d/2) is the surface area of a unit d-dimensional sphere. Upon elementary transformations and exploiting 
the integral representation of the hypergeometric function. 



F{a,b, c; z) 



r(c) 



r(6)r(c-6)7. 



dtt''-\i-ty-''-\i-tz)-'', 



it is possible to obtain 



d(t)I[ct>,p] 



r(d/2) 
r(f - 1) 



turn 

turn 
pq 



F 



F 



1 d 
2'2' 



1,2- 



/pgy 

\uin ) 



d 3 / Lom 



2 2' \ pq 



{c>b>0), 

(ujm > pq) 
{pq > com) 
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(B6) 



(B7) 



and 



dffl sm (h — — — - — -4 — . 

r('^/2) 



(B8) 



If we assume, as is the case for a degenerate electron gas, that ksT is much lower than the Fermi energy then the momentum 
integral can be simplified with the substitution -^ — — (5(e — ep)- The final expressions, at this low temperature, are then, 



Rexn(q,t^) = 



7r'*/2r(d/2) 
7n2^-''p^-"+" 



1-F 



1 i dJppqV 
' 2' 2'\ujmJ 



2r(d/2) [cum 
7rrf/2r(d/2) ] r(-f + d/2) [ppq 



d 3 j uom 
^'^"2'2'1^^ 



(a;m > ppq) 
[ppq > ujm). 



(B9) 



The static limit of ReXn is thus identified with the function, _B„(q), for a highly degenerate electron gas. 

APPENDIX C: GENERALISED SOLUTION FOR THERMODYNAMIC FLUCTUATIONS 

The macroscopic fluctuations of Sn, SJ and SE, as defined by Eq. ([12|), can be obtained by taking the relevant momentum 
moments of the microscopic equation of motion which we rewrite here in terms of q. 



(^-Ae){p+iq|5>!p-iq) + A/(/(q) = -^ f (p + iq|^'^p|p_ iq) + ^ 



«„ + p.M<,) + (^ + |j-.)^ 



(Cl) 



However, since we assert that the macroscopic quantities are conserved before and after each collision, we must obtain 
two sets of expressions for each term. The first set, determined by the dynamics after each collision process (i.e. ri ght- hand 
side of Eq. (|Cl|)) is given by Eq. ([l2t). The second set follows from solving the entire microscopic equation, Eq. (Cl), for 
{p+|qj(5p|p— |q) and then taking moments, i.e. 




P[Ae] 
P[pAe] 

L2r7T -I 



/ 



PMq) + p-<5v(q) + (^ + i;r-M) 



STid) 



P[p(5/i(q) + p . 5v(q) + (1^ + i;^ - /i) 



STid) 



\ 



\p[£-(S^^(q)+p■5^.icl) + C- + ^^-^^) 



5T(q) 



(C2) 



T 



where 



^5] = 



(2^) 



,d A/ 



Ae d) - Ae ■ 



Upon expansion of partial fractions we can obtain 



Co 




Bo + Co Co-^ 



-B2+C2 

2m 



+ (£-m)(5o + Cu) 



-t/(q) ^iCo + -^ 
\ C2_ li^r 

V 2m / 



Coti 



-B2 + C2 

2m 



Gof; i^<i+(£-M)Coi 



£2. 

1m q^ 



1m ^ 8m 

B4+C4 I iqL _ „'| B2+C2 
4m^ ' V 8m i^l 2m 




(C3) 



(C4) 



where we have, for convenience, taken the component of momentum density along q. In addition we have assumed that the 
direction of the momentum density lies entirely in the q direction, a property that fails for Non-Newtonian fluids. With this 
expression, together with Eq. (|l2), it is possible to eliminate the conserved quantities to solve for the thermodynamic variations. 



Do 



Co 



m) A) 



IT iTo-^ zmiT ^ am ' ■' it 



£2. 4. (j_ 

2m ' \ 8m 
C2_ 



D2 
2m 



C2 



Imq^ IT 



Am-i ~^ V 8rr 



^^)Tr. 




(C5) 
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FIG. 1: Plot of Re(e) against u/Ep for a three-dimensional system. The parameters are r^ = 3, g = (Q.2)kp , t — 16/-EJ^' 
andT = OMEp/kB. 







FIG. 2: Plot of Im(e) against lo/Ef for a three-dimensional system. The parameters are the same as in Fig. hi 
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FIG. 3: Plot of Re(e) against lu/Ef for a two-dimensional system. The parameters are Ts — 3, q — {0.2)kp , r = 10/ Ep 
and T = 0.01EF/kB. 
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FIG. 4: Plot of l7Ti(e) against u/ Ep for a two-dimensional system. The parameters are the same as in Fig. H. 
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FIG. 5: Three-dimensional plasmon dispersion. The parameters are rs = 3, r = ?>/Ep and T = O.Ol-EF/fcs- The arrow 
indicates the plasmon frequency, Wp = 47re^no/m. 
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FIG. 6: Two-dimensional plasmon dispersion. The parameters are r^ = 3, r = 5/Ep and T = O.Ol-EF/fcs 



